Engineering interband transport by time-dependent disorder 
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We show how the evolution of atoms in a tilted lattice can be changed and controlled by phase 
noise on the lattice. Dependent on the characteristic parameters of the noise, the interband transport 
can either be suppressed or enhanced, which is of interest for high precision control in experimental 
realization with Bose-Einstein condensates. The effect of the noise on the survival probability in 
the ground band is summarized in a scaling plot stressing the universality of our results. 
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Ultracold atoms in optical lattices have opened a pos- 
sibility to investigate effects which were not previously 
observable on ordinary matter crystals, at least not in 
such a controlled way, like Bloch oscillations, Wannier- 
Stark ladders, and Landau-Zener tunneling [l|-l7|. Load- 
ing Bose-Einstein condensates (EEC) into optical lat- 
tices provides an optimal control of the studied systems 
by optical means. More complicated potentials can be 
realized by adding further lattice beams [1, . By super- 
imposing laser beams from different directions and with 
different frequencies, it is possible to generate various lat- 
tice geometries [1, [9l-[lll. \W{ . The question arises how to 
control the dynamics of particles by quasi-periodic poten- 
tials (possibly time-dependent or even stochastic ones). 
From the theoretical point of view, a variety of phenom- 
ena is expected to occur in these systems, such as An- 
derson localization [ill [T^j and the quantum transition 
to the Bose glass phase originating from the interplay of 
interaction and spatial, but static disorder 
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In this Rapid Communication, we investigate the 
Landau-Zener tunneling of a BEC from the ground band 
of a lattice into higher bands and, finally, into the ener- 
getic continuum. We show that control on the interband 
coupling is not only achieved by varying the lattice pa- 
rameters, the initial conditions 0,01, or by stimulated 
transitions [lol |. but also by additional noise which pro- 
duces a time-dependent disordered potential [lB|- More 
specifically, we study the time evolution of the survival 
probability of a BEC loaded into a quasi-lD geometry 
built of a stochastic potential and a static Stark force. 
While preliminary results reported in [l6| imply that ex- 
ponentially correlated noise is not efficient in controlling 
the interband decay, we present here a universal predic- 
tion based on harmonic noise with a well-defined fre- 
quency range. By adapting the noise frequency to the 
system's scales we demonstrate that stochastic noise can 
be used to engineer the decay very efficiently. 

We focus on the dynamics of a sufficiently dilute non- 
interacting BEC. In the absence of atom-atom inter- 
actions, the temporal evolution of a BEC in a tilted 
one-dimensional lattice is described by the single-particle 
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where M is the atomic mass and F the Stark force. 
Vs (x, t) is a stochastic potential given by the following 
time-dependent bichromatic lattice 



V,{x,t) = aV 
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It consists of two spatially periodic potentials with in- 
commensurate spacings, c^l for the "reference" lattice 
and c?L = c^l(%/5 — l)/2 for the secondary lattice. The 
time-dependence arises from the time-dependent stochas- 
tic phase (^{t). For constant 0, this is the system with 
a quasi-disordered lattice which was realized by Roati et 
al. [ll[. The necessary phase control is easily achieved 
in experiments for the time scales studied here (about 
10^s...0.1s) using modulated lattice beams [l^l- The 
depths of the two lattices are considered to be compara- 
ble and for convenience equal amplitudes aV are chosen, 
where V is the lattice depth of the reference lattice and a 
is a renormalization factor (which will be defined below) . 
The recoil energy E^^c = Prec/2A/, with p^^c = 7r?i/dL, 
is the characteristic energy scale of the optical reference 
lattice 0, and Vq = V/E^cc and Fq = Fdi,/E^cc are di- 
mensionless quantities in this energy unit. We assume 
a typical experimental situation [J [j, @] , for which the 
initial state is a condensate's wave function relaxed in a 
harmonic trap and then loaded adiabatically into the lat- 
tice given by Vs{x,t = 0). The trap is shallow such that 
many lattice sites are initially populated. It is turned off 
for the further time evolution. 

For the original Wannier-Stark problem including just 
the first term on the r.h.s. of Eq. ([2]), the wave packet 
prepared in the ground band oscillates periodically in 
time with the Bloch period T-q = 27Th/ (FqE^cc) and fre- 
quency LdB = Stt/Tb [17]- However, for sufficiently large 
F, the BEC undergoes Landau-Zener tunneling after 
each Bloch period at the band edge, where the band gap 
has its smallest value [1, [l^ . Such a phenomenon gives 
rise to a step-like structure in the time-resolved survival 
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probability Psurit) of the atoms in the ground band 0. 
Following 

0, S, [Hi, we calculate Psur(<) in momentum 
space by projecting onto the support of the initial state: 
Psurit) = (ipx|V'(Px, OP- The upper bound can be 
safely extended to infinity since the tunneled wa ve p acket 
moves into the direction of negative momenta [18|- The 
lower bound is given by the parameter pc, which can be 
adapted to stabilize the results with respect to temporal 
fluctuations. While the true decay starts after a traver- 
sal of half the Brillouin zone (the Brillouin zone has a 
width of 2pi.oc), we choose either pc = 3prcc or pc ~ 5prcc, 
which effectively means that we measure the decay with 
a delay of Tb or 2T^ respectively. As long as the mean 
momentum evolves linearly in time with a variance sig- 
nificantly smaller than the width of the Brillouin zone, 
we can resort to the acceleration theorem [l[ to relate 
both ways of measuring. In the following, we introduce 
the type of stochastic process we use to generate (/)(t) and 
show universal properties of the survival probability as a 
function of rescaled parameters. 

In general, colored noise can be defined by a suitable 
spectral distribution. The simplest type can be produced 
by a harmonic oscillator of angular frequency wq and 
damping rate F, driven by Gaussian white noise. The 
result is harmonic noise [l9|, which is represented by a 
two dimensional Gauss-Markov process taking the form 
of the following stochastic differential equations for 4>(t): 

(j) ^ V (3a) 
V ^ -2Tv - ujIcP + S,{t) . (3b) 

^(i) is Gaussian white noise with (^(t)) — 0, (^(i)^(s)) = 
ATT5{t — s), and T measures the strength of the noise. 
F and ojq have the units of Tg ^ and ojb respectively. 
For notational reasons, we express both of them in the 
unit of ^- The variance of the noise is given by 
((/)^(i)) = T/wg, where (•) is a time average over a suffi- 
ciently large interval (in our simulations ~ 20 Tb). The 
spectral distribution for the harmonic noise is 

2Fr 

^ ^'^^ " 7r(4F2cj2 + (cj2 _ ^2)2) ' (4) 

This distribution incorporates two important regimes: 
(i) slowly-varying noise when the oscillation frequency is 
much less than the damping rate, i.e., 2r^ ^ luq, and (ii) 
fast noise, when 2F^ ^ ujq. The spectral distributions 
corresponding to these limits are shown in the insets of 
Fig. [TJ In the regime (i), the noise recovers the expo- 
nentially correlated noise [l^, [2^ with a Lorentzian-like 
distribution peaked around zero frequency (see Fig.[T][a)). 
The spectral distribution in the regime of fast noise (see 
Fig. [TJb)) peaks at ujq = -^/cjq — 2F2 with a finite width 

AujQ ~ 2'\jTy^Q — F^. To control the evolution in the 
potential of Eq. ([5]) the most relevant noise realization is 
this latter regime (ii) , for which the above mentioned en- 
ergy scales can be matched with the scales of the original 
Wannier-Stark system. 



For very fast noise, the particle effectively averages 
over the time-dependent potential. Assuming a Gaus- 
sian distribution of the fast varying phase with variance 
{(j)^) and integrating over the phase as random variable, 
an effective potential is calculated as 



Kff (x) = 




where /3 = exp {-2{(j)'^)) . The effect of fast noise is hence 
to renormalize the amplitude of the second lattice by a 
factor (3. The parameter a is introduced to compare 
better the dynamics in the potential given by Eq. ([2]) 
with the dynamics in the reference system, in which just 
V (x) = V siia^ (nx/di^) is present. For this, a is chosen 
such that the spatial second moments of V{x) and Vcb{x) 
are equal [3 ■ 

We study the temporal behavior of the survival proba- 
bility Psurit), which is the main experimental observable 
for our system 0, @, [13 > in a- broad range of noise pa- 
rameters. We use a reference system with typical exper- 
imental values of Vb = 2.5 (giving an average band gap 
AE 2.5i?iGc) and _Fo = 1-5 [l,[ll. wq and F, are chosen 
in the range of O.OI/Tb and 300/rB, which covers both 
regimes of fast and slow noise. The stochastic nature 
makes it necessary to average the results over a sufficient 
number of statistical realizations. For the parameters in- 
vestigated here, 20 realizations for fixed noise parameters 
turned out to well stabilize statistical fluctuations. The 
lattice parameters Vq and Fq can be adapted to fulfill the 
condition of resonantly enhanced tunneling (RET) in the 
reference system, i.e., Fq « nAE / Ej-^c, with an integer 
n [1, [l3, [3 , or away from this special condition (non- 
RET) . Since at RET the transport from the ground band 
to higher bands is already enhanced by energetic quasi 
degeneracies [l3l, and the noise is likely to drive the sys- 
tem out of RET conditions, the decay is expected to de- 
grade for typical parameters in this case [la]- For our 
non-RET conditions Vq = 2.5 and Fq = 1.5, it is shown 
in the following that a faster decay can be easily induced 
by harmonic noise, since there is no competition between 
the two effects (enhancement by RET and noise). 

The stochastic potential Vs {x, t) is determined by the 
spectral distributions S{uj) of the noise, which are de- 
picted in the insets of Fig.[TJa) and (b) for the two differ- 
ent regimes of noise, respectively. For a small oscillation 
frequency loq = 0.1/Tb compared with the damping rate 
F = 5/Tb, S{uj) has a very narrow peak at zero frequency. 
The phase itself is slowly varying with time in case (a) , 
whilst it shows faster fluctuations in case (b). The effect 
of the noise on the temporal evolution of Psur('^) is com- 
pared in Figs. (Ha) and (b) to the one for the noise-free 
reference system with its characteristic step- like structure 
(sohd lines) d [13. As seen, the harmonic noise tends 
to wash out the step structure after a few Bloch oscilla- 
tions. Moreover, it leads to a systematic enhancement of 
the tunneling rate for the largest noise amplitude T in 
both cases (a,b). In the regime of fast noise (Fig. [ijb)) 
the tunneling rate is always enhanced with respect to the 
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FIG. 1: (color online) Psur(t) Ki = 2.5, Fq = 1.5 in the 
presence of harmonic noise with F = 5/Tb, and ojo = O.I/Tb 
((a) - non-universal regime, c.f. Fig. [2]at low values on the x 
axis) or wo = IO/Tb (b) and for T = 0.01/T| (green dashed 
lines), T = l/r| (blue dot-dashed lines), T = IO/tI (red dot- 
dot-dashed lines), and T = 100/r| (dotted lines). The result 
for the reference system is shown by the thick solid lines. The 
corresponding noise spectra S{iu) are given in the insets for 
the case of T = 10/r|. 



reference system, and here it is essentially independent 
of T (sec also Fig. [3] below). 

We find that, for a broad range of noise parameters, 
our main observable, the survival probability Psur(^o) a-t 
fixed time to, obeys a scaling relation as a function of all 
parameters. Our results are similar for other choices of 
^0 ^ ^Tb; and we stick in the following to to w 6Tb, a 
typical experimental observation time reported, e.g., in 
Q . In Fig. [21 the noise frequency wq is varied by keeping 
the variance of the noise fixed. Rescaling now wq by ujb, 
by the standard deviation of the noise y/ (0^) and by the 
ratio of the two lattice periods d-^/di^, all the results of 
Fig. [2] show a similar behavior with a dip at ujo\/ (0^) ~ 
ijjBdi^/(2di^). First this value of energy lies just above 
the average band gap, i.e., wqa/ (0^) ^ AE « 2.5i?i.cc, 
which turned out to be a necessary condition to observe 
the dip. Secondly, it implies that the optimal time scale 
of the noise is directly related to Tb for observing an 
enhancement of the decay 

The additional factor d^/di^ takes account of the dif- 
ferent periods of the two lattices, and the scaling was 
checked for various incommensurable ratios dj^/d]^ (not 
shown) . Faster noise is less efficient again since it leads to 
the regime where the effective model applies, see Eq. ([5]). 
Here the noise effectively averages during the evolution. 
Indeed at very large wq, -Psur(^o) lies on top of the corre- 
sponding value (shown by horizontal lines) calculated us- 
ing the effective potential. The observed scaling relation 
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FIG. 2: (color online) Psur(io ~ 6 Tb) for Vo = 2.5, Fo = 1.5, 
F = 5/Tb, for fixed noise variances (0^) — 0.25 (correspond- 
ing to a ~ 0.85, open circles), 1 (a « 0.99, filled squares), and 
10 {a ~ 1, open squares). For comparison we show Paur(to) 
for the respective effective models by the dash-dash-dotted, 
dot-dot-dashed, and dashed lines (in the same order and color 
of the symbols), for the reference system (thick solid line), 
and calculated using a deterministically oscillating phase with 
Al/2 = T/uja = 0.25 (pluses), 1 (stars) and 10 (crosses). 



around the dip, as our main result, drastically reduces 
the parameter dependence of the system. It helps pre- 
dict and control the behavior of the survival probability 
by a simple choice of right combinations of the a priori 
many parameters (Fq, Vq, Ti 2^)- 

Deviations from the scaling are expected when the 
noise frequency wq becomes comparable with its spectral 
width, which is determined by F. This implies a crossover 
from the spectral distribution shown in Fig. [Ijb) to the 
case of Fig. [IJa) . Then the noisy phase changes slowly 
in time, and for T < 10/r| cannot assist the tunneling 

anymore, as observed in Fig. [ija). Here uJoy^ (0^) ^ wb 
and the interband tunneling depends very much on the 
strength T of the noise, cf. Fig. [51 since single noise 
events can dominate the evolution. To better understand 
the regime where the scaling breaks down, it is helpful 
to study more explicitly the dependence on {(f)'^), which 
is done in Fig. [21 We distinguish two cases: in (a) the 
data is produced as previously described by adapting a 
to compare with the reference system (a is calculated as 
a function of /3 = exp(— 2(0^))); in (b) we keep a fixed to 
understand the effect of the absolute value of the lattice 
depth (loosing, of course, the meaning of a reference sys- 
tem). The data collapses onto a single curve described 
by the effective fast noise model of Eq. ^ (dashed lines 
in Fig. [3]) for large frequencies uJq/T > 100. The behav- 
ior of the various curves in (a) is similar (approximate 
scaling) for 1 < wg/T < 100, provided that T < 1/T|. 
As shown also in Fig. [IJa), for a large T > lOO/Tg even 
the slowly varying noise is strong enough to enhance the 
decay below the value predicted for the static reference 
system. For the latter case, the scaling relation obviously 
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FIG. 3: (color online) Psuv{to ~ 6 Tb) for = 2.5, Fo = 1.5 
vs. {4>^)~'^ = oJo/T: (a) adapted a / 1; (b) fixed a = 1. 
Other parameters: T = 5/^8, and T = 0.001/Tji (filled cir- 
cles), T = O.OI/tI (open circles), T = 0.1/r| (filled squares), 
T = l/rl (open squares), T = 10/r| (filled triangles) and 
T — lOO/Tg (open triangles). Psui (to) of the corresponding 
static reference system is shown by the solid line in (a) and 
for the effective model with /3 = 1 by the dashed lines in (a,b). 



is bound to fail. In Fig.[3jb), where we do not correct for 
the change in the relative height of the lattice amplitudes, 
the just described trends in (a) are shifted and the scal- 
ing is harder to appreciate (except for the very fast noise 
again, where our data always follows the expectation of 
the effective static model). 

As a final benchmark, we compare to a periodically 
oscillatory phase (j){t) = cos{ujdt + if), with being 
its amplitude and Wd its oscillation frequency. To model 
random noise fluctuations the data are averaged over the 
parameter ip (picked randomly from a flat distribution). 
The spectral distribution of this function is given by 



S'd (w) = ^ {(5(a; — Wd) + (5(cj + cjd)}. Practically, the 
delta functions gain a finite width Aaj,„i„ = "^here 
Ttot is the total numerical integration time, fke spec- 
tral distribution of the harmonic noise in the fast noise 
regime, which very much peaks at wqj is then compara- 
ble to the delta-like spectral distribution around Wd of 
the oscillating noise. Therefore, for Wd ~ '^o ^ T and 
Ad = \J2T Ilo\, the periodic phase (averaged over ip) has 
a similar effect as harmonic noise. Examples of results 
are shown in Fig. [5] for A^/2 = 0.25, 1 and 10, respec- 
tively. As expected, for large the symbols lie again 
on top of the prediction of the effective model. 

We investigated the impact of noise on the Landau- 
Zener tunneling in a one-dimensional Wannier-Stark sys- 
tem. By a proper scaling of the data with parameters, 
we can universally characterize the effect of the noise. 
This central result is summarized in Fig. [2] Our predic- 
tions show that time-dependent noise in a bichromatic 
lattice provides a further handle to control the transport 
from the ground band to higher energy bands. We de- 
scribe a first step in the direction to push investigations 
of the static regime into the realm of dynamical disor- 
der [2l[ . A natural extension would be to study the case 
of the simultaneous presence of noise and atom-atom in- 
teractions in our Wannier-Stark system, complementary 
to expansion experiments in flat non-tilted potentials in 
optics [H and with BEG [H]. 
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